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We revisit the issue on signatures of pre-inflationary background anisotropy by considering the 
quantization of a massless and minimally coupled scalar field in an axially symmetric Kasner back- 
ground, mimicking cosmological perturbations. We show that the power spectrum of the scalar field 
fluctuation has a negligible difference from the standard inflation in the non-planar directions, but 
it has a sharp peak around the symmetry plane. For the non-planar high-momentum modes, we 
use the WKB approximation for the first period and the asymptotic approximation based on the de 
Sitter solution for the next period. At the boundary, two mode functions have the same accuracy 
with error of O(Hifk). We calculate the approximation up to the order of (Hi/k) 6 and show that 
r^j . the power spectrum of the scalar field fails to get corrections until we execute the approximation up 
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I. INTRODUCTION 



Inflation has become one of the paradigms of modern cosmology. First of all, inflation elegantly solves many 
I , problems which are present in the standard Big-Bang model such as the horizon and flatness problems. Second, it 
jHi accounts for the origin of the large scale structure of the universe in terms of the quantum fluctuations originating 
from the adiabatic vacuum structure in early universe. According to the cosmic no-hair conjecture, the exponential 
expansion of the universe during the inflationary period erases any classical memory Remarkably, the nature 
of the primordial fluctuations is understood in terms of symmetries of the de Sitter spacetime 0. In general, we 
^ , need n-point correlation functions to characterize the statistical nature of primordial fluctuations. However, these 
symmetries lead the power spectrum of a scale invariant form: 



p(0) _ m 
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where Hi is the Hubble parameter during inflation. These predictions from symmetries are robust and universal 
in inflationary scenarios. In fact, the above predictions have been confirmed, e.g., by the measurements of cosmic 
microwave background (CMB) 0. 

As the observational precision increases, we have to go beyond the power spectrum to look at fine structure of 
the primordial fluctuations. Since in the realistic inflationary universe the symmetries of the de Sitter spacetime do 
not hold exactly, violation of them provides a measurable effect. The violation of the temporal de Sitter symmetry 
leads to a slight tilt of the power spectrum, which is characterized by a slow roll parameter. Similarly, the violation 
of the shift symmetry leads to non-Gaussianities Q, for which we need to prove higher correlation functions. It is 
known that observable non-Gaussianities do not appear in the simplest single field inflaton with a canonical action. 
Substantial progresses [f| have been made in understanding the enhancement in non-Gaussianity for several variants 
of the standard scenario, involving having multiple scalar fields, non-canonical action for the scalar field, introducing 
higher derivative terms in the action, or having a non-standard vacuum state. 

One may naturally expect violation of the spatial de Sitter symmetry, which would lead to the statistical inho- 
mogeneity 0. Another possibility is introducing spatially homogeneous models violating the spatial isotropy (7rll0j. 
where the Copernican principle is kept since there is no privileged positions in the universe. Instead, in this model, 
the universe has a privileged direction. To reconcile with the spatial homogeneity, the axis should be the same all over 
the universe. The existence of a privileged direction which is global over the universe seems unlikely to happen. It 
may be possible only when the speciality comes from some quantum mechanical process in the very early universe. It 
is known that a primordial anisotropy in the metric will be washed away during inflation, hence it is interesting to find 
remaining signatures of the strong anisotropy in the early spacetime metric in the perturbation spectrum. From the 
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observational point of view, a lot of anomalies indicating the statistical anisotropy are reported although its statistical 
significance is uncertain. Unfortunately, it turns out that many of these models suffer from either instabilities, or a 
fine tuning problem, or a naturalness problem (Tlj . 

We start from the discussion on the evolution of anisotropic universe in the Einstein gravity minimally coupled 
to a massive scalar field, where the scalar field plays the role of the inflaton. To obtain a sufficiently long period of 
inflation, one usually imposes the slow rolling condition, 4>o 3> Mp, where 4>o is the initial value of the inflaton. Under 
the assumption, as discussed in Ref. (l2j . the background metric can be approximated by the Kasner spacetime with a 
positive cosmological constant, A(= 3iJ; 2 ). Among all, we are mainly interested in the the regular Kasner spacctimcs 
with two dimensional axial symmetry with metric, 
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The spacetime has a privileged axis £3 orthogonal to the symmetry plane. The spacetime has a Rindler-like event 
horizon at r = since gn(= 322) approaches to a finite value and 533 goes to zero linearly. The region r > 
is connected to the region beyond the horizon where the coordinate r plays the role of a space coordinate and a 
timelike singularity is located. As shown in Ref. [l3j], the geometry is closely connected with the anti-de Sitter black 
brane solution discovered by Lemos [Hj]. There are various reasons why this spacetime is a good testing ground in 
analyzing the properties of anisotropic universes. First, it bears various important features of the whole anisotropic 
universes including large anisotropy at r = 0. Second, it is known that this anisotropic inflation is an attractor Q 
in supcrgravity with a wide range of gauge kinetic functions. This fact will be related to the well-known instability 
of gravitational wave modes in Kasner universes [tj [l(| , since the instability becomes controllable only for the above 
case ([2]). Third, only in this symmetric case of all aniso trop ic expansions, we can impose proper anisotropic vacuum 
state in terms of the zeroth order WKB approximation [15[ . 

If one uses the Sasaki-Mukhanov variable, except for the complication due to the mixing of tensor-scalar modes, 
the evolutions of the metric perturbations are not much different from that of a massless scalar field [l|| [l8|. In 
Ref. (l2l ] , it was also shown that the mode mixing modifies the power spectrum only by a proportionality factor with 
a small correction term of order ml Hi, where m is the mass of the inflaton. Therefore, as a formulation level, it is 
good to deal with the scalar field rather than the metric perturbation itself. Hence in this work, we are interested 
in the evolution of a massless, minimally coupled scalar field propagating on the background anisotropic universe ([2]) 
with action 
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This scalar field is not the inflaton mentioned above but just a mimic of metric perturbations. Hence it is assumed 
that its size is small enough so that it does not modify the background dynamics. We focus on the quantization 
procedure and the evolution in the anisotropic stage. The effects of a scalar field onto the geometry of the anisotropic 
Universe were discussed in [19j . The canonical quantization of the scalar field can be done in the standard manner. 
The vacuum is chosen at the initial anisotropic era: r — > +0. Usually, in a time-dependent quantum harmonic 
oscillator, an adiabatic process implies the one where the potential changes slowly enough compared to its size, and 
the time evolution can be obtained from the WKB approximation. In the standard inflationary models, an adiabatic 
vacuum in the early universe is also defined in the same way. However, in anisotropic space-times, the situation gets 
worse for two reasons. First, most modes do not take oscillatory form, therefore, the adiabatic vacuum does not exist. 
Second, there exists the instability problem for the tensor modes of gravitational perturbations. Only in the specific 
solution (|2l . the two problems disappear and the adiabatic vacuum exists, which was called an anisotropic vacuum 
in Ref. [la ]. It is important to note that, in this scheme, the entire information of anisotropy (and therefore, the 
cosmological structure at the present universe) is encoded in the vacuum state of the theory similarly to the case of 
the Robertson- Walker space-time. 

In general, the evolution of the scalar field in this spacetime cannot be solved exactly. To have an analytic form of 
the spectrum, we develop an appropriate approximation method. Usually, the WKB approximation and the de Sitter 
solutions are commonly used at the early and later stages of evolutions (l2l . [TEl [l6j . For ordinary high momentum 
modes with k 3> Hi, the WKB approximation is valid. However, at the very early stage of the anisotropic expansion, 
the validity condition for the WKB approximation changes to Hi, where k% is the comoving momentum along 

the special axis. Therefore, there exist non-trivial high momentum modes (the planar modes, k% ~ Hi and k Hi) 
which cannot be dealt with the WKB approximation but need separate treatments. 

Docs the present approximation successfully predict the power spectrum for the present universe? For the observed 
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wave length scales smaller than 10 4 Mpc, the mode with comoving momentum fc Q b s can be seen if it satisfies [HI, HH 

k , >k , = H e N ~^ ( Tr \ ( 1Ql6GcV ^l 2 (4) 

where fc C rit, -^V, Tr, and Vl n f are the critical observable comoving wavenumber, the number of e-folding, the reheating 
temperature and the value of inflaton potential, respectively. Thus, for inflation in the grand unified theory (GUT) 
if e-folding is around N > 64, the modes may satisfy fc Q b s ^ Hi. Therefore, most observable modes in our universe 
marginally satisfy the WKB approximation condition. In this sense, the present approximation can successfully expect 
the evolutions of all non-planar modes. For a sufficient long duration of inflation, say N 3> 64, all observable modes 
satisfy fc 3> H and the WKB approximation holds well. However, in this case, it will be very hard to find any signature 
originating from the initial anisotropy since most correction terms will be suppressed by some power of Hi/k ~ e 64_JV 
except for the planar modes. If N ~ 60, the modes which belong to the wavelength scale between I0 2 Mpc and 
I0 4 Mpc, will bear much information about the anisotropic prcinflationary memory, which must be already detected. 
This may constrain the e-folding of the inflationary expansion from being smaller than 64. Assuming N > 64, it is 
not difficult to satisfy the WKB condition for the range of observable scales. However, the visibility of the anisotropy 
at these scales is not guaranteed. 

In this work, we show that the relic of the prcinflationary anisotropy appears to be very difficult to be observed in the 
non-planar modes except for the special case that the residue of the comoving momentum satisfying Hi < k < lOiJ; 
are located in the observable range. This is simply because the non- vanishing correction to the power spectrum comes 
from very high order perturbative corrections and therefore, its size is suppressed by high power of Hi/k. Pcrturbative 
calculations on the power spectrum were done previously in Ref . [l2l . [Til [iH ] , where the authors are missing doing 
perturbative approximation for the asymptotic solution around the de Sitter one. Still, there are possibilities that the 
anisotropic effect can be detected only from the power spectrum for planar modes. 

In Sec. II, we quantize the scalar field and gives the initial condition for the scalar evolution by setting the mode 
solution to be purely positive frequency at early times. In Sec. Ill, we introduce the approximation method used in 
this work and show the resulting power spectrum of the scalar field up to 6 th order of the expansion parameter. We 
show that the non-trivial correction to the power spectrum from initial anisotropy survive only from that order. In Sec. 
IV, we present the power spectrum for the case of a planar modes where the contribution from the initial anisotropic 
expansion of the background geometry appears in its zeroth order approximation. In Sec. V, we summarize the 
results. We attach a few Appendices where the details of the calculations are given. 



II. QUANTIZATION OF THE SCALAR FIELD ON ANISOTROPIC VACUUM 

The canonical quantization of the minimally coupled massless scalar field with the action ([3]) in the anisotropic 
spacetimc ^ is done in the standard manner: 

4> = j d 3 k(^u k a k + u^al^j, (5) 

where the creation and annihilation operators satisfy the commutation relations [aid,^ ] = <5(ki — k2) (others are 
zero) and «k = e lkx </>k/(27r) 3 / 2 . The details of the quantization process depend on the choice of the mode function 
Uk- We normalize the mode function as 

Writ- = ± (6) 

For the later convenience, we introduce a dimcnsionless time x by 

Sinh ^-sim^= e ^' < 7 > 

where e denotes a small expansion parameter which will be specified later. The arrow of time for x is inverted since 
it varies from oo to 0+ as the comoving time r increases from 0+ to oo. The asymptotic limit with r 3> H^ 1 is now 
equivalent to ex -C 1. Therefore, the asymptotic expansion in the long time limit corresponds to the e expansion. 
The equation of motion for the scalar field is written in the form of a time-dependent oscillator 

(^2 +«kW 2 )0k = O, (8) 
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where the dimcnsionless frequency squared is 



kW V3^J (l- e -2-)4/3 " g ^sinhexy' ^e 2 - - 1 + J ' l j 

Here we define a scaled wave-number and a measure of planarity of a given mode by 

i-*»±, r-|, (10) 

where fc 2 := A: 2 + fc| + fc| = fc 2 ^ + fcf. Later in this paper, we omit k in the frequency squared fi 2 . for simplicity. In 
terms of x, the normalization condition ([5]) changes to 

if 

^kd^-(^ k )^ = - — . (11) 

Here, we use dr/dx = —e 3a x e/(3Hi). 
The power spectrum is defined by 

(O|0 2 |O) := J d\nk J P = |^ k | 2 . (12) 

In contrast to the case of the standard de Sitter universe the direction dependence would be included in the power 
spectrum. The vacuum is chosen at the initial anisotropic era: r — > +0 to satisfy flk|0) = 0. For this purpose, we 
choose the solution to be purely positive frequency mode with respect to r at the early stage: 

hm k = 1 e ^+ # (13) 

^° a/2 x 2 2 /3fc 3 



where we follow the normalization ([5]) and 
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III. NON-PLANAR HIGH-MOMENTUM MODES 



As mentioned in the introduction, we use the WKB solution at early times (x > £*) and asymptotic solution at 
later times (x < £*). At the matching time x = x*, the accuracies of the two are equal. Later in this paper, we set 
= 1 by choosing e appropriately. We also assume that k is larger than one, which will be satisfied with the modes 
we are interested in at the present approximation. 



A. WKB solution at the early times 



The WKB approximation is valid up to the order of the correction term if 



e(x) 



dn 2 {x) 

dx 



n 3 (x) 
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I V i 



1/3 



-2ex 



1 . 2 

T, tanh ex 



vT 



-2sx 



1 



3/2 



< 1, 



(14) 



where e plays the role of the adiabatic parameter. For most modes, the adiabaticity condition is valid if ex 3> 1 
because of the exponentially decaying factor e~ 2ex . For a small ex, the condition is also valid for non-planar high 
momentum modes if x 3> k~ 3 . In the asymptotic regime with < x < k~ 3 , the adiabaticity condition fails to be 
satisfied. On the other hand, the asymptotic (small ex) approximation is valid if x <C e^ 1 . The asymptotic expansion 
in this case will be dealt in the next subsection. 

The WKB solution satisfying the initial condition f| 13[) is given by 



'WKB 



i 



2Q(x) 



exp 



i / dx'Cl(x') + it/i 



(15) 
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where we choose the positive (negative) frequency mode with respect to t (x) and ip is a phase factor. Q, satisfies the 
nonlinear equation 

n(x) 2 =n(x) 2 c(x) 2 ; c(x) 2 := 1 - iW + -( V (16) 

The WKB wavcfunction is expanded up to an enough adiabatic order, to validate our matching scheme with the 
solutions in the asymptotic region. Noting Eq. (|16[) . the next order of the approximation will improve the accuracy 
by the order 

Ewkb{x) ~ <x) 2 ~ (17) 



1/3 

fc 2 £C 2 /3 T fc4 a ,4/3 "1 3 p" 
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For example, to the order fc 4 , we get c(x) 2 = 1 — fr2 2 2 / 3 + fr4 ^ 4/3 + 3r 3 1 1 — • , where the last term has the 

mixed order of magnitude e/k 



B. Later time solution 



In the later time limit, the space-time approaches that of the de Sitter spacetime. Therefore, it is natural that the 
zcroth order solution is just the well-known scalar field mode solution in the de Sitter spacetime. The higher order 
corrections are based on the asymptotic approximation based on the limit a<Cl. By series expanding the frequency 
squared, we approximate the equation of motion (j5J) to be 

(d 2 °° \ 

(^+£^U = °. as) 

\ n=0 / 

where the order by order corrections of the frequency squared are 

P 2k 2 q 2 _ 2fcW/ 3 _ 8Px 5 / 3 

V °~9^' ~~ 27 ' V3 -^29~ r " (19) 

with q 2 = 1/3 — r 2 . For higher orders than these, refer Eq. (|B2[) in Appendix [Bj Note that each Vi is of the same 
order as Vq in this expansion. This is why we prefer x to r. We present the corrections up to the 6 th order in e in 
Appendix [B] This is because nontrivial corrections to the power spectrum of the isotropic case exist from this order. 
The general solution to the differential equation (|18p is given by 

cf> = A + u{x) + A_v{x) , (20) 

where u and v are the positive and negative frequency modes, respectively. We get the approximate solution by series 
expanding the modes u — uq + eut + e 2 U2 + ■ ■ ■ and v — vq + ev\ + e 2 v 2 + ■ • • and then solving the equation of motions 
order by order. The zeroth order equation satisfies 



2 



0. (21) 



dx 2 9a; 4 / 3 

Its solutions are nothing but the de Sitter mode solutions 

u (x) := ( - 1 + ikx 1 ^^ 3 , v (x) :=(-!- ikx 1/3 ^j e "^ 1/3 = u *(x), (22) 

where the Wronskian becomes uq(x)v' (x) — u' (x)vq(x) = — 2ik 3 /3. This Wronskian condition will be kept even if we 
get higher order corrections. 

The equation of motion for the first order correction term is 

^ + ^)^ = -Wo(4 (23) 

The left hand side of the equation of motion is the same as that of the zeroth order and — V\{x)4>q{x) plays the role 
of a source. Its solution is formally given by the integral, 

M*) = - I " G(x,y)V 1 (y)My)dy, (24) 
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where the Green's function is, 



G(x,y) 



2ik 3 



[u (x)v (y)9(y - x) + u (y)v (x)9(x - y)] 



(25) 



One may choose the upper bound of the integration x c to be any value greater than x formally. In fact, the choice of 
x c changes only the homogeneous part of the solution (explicitly 4>o), which is not relevant in doing the perturbative 
higher order calculation. Therefore, we simply drop all x c dependent term after the integration so that the first order 
corrections ui is the complex conjugate to v\. Inserting 4>o(y) — > Uo(y) or — > v (y) we get 



i>*(x) = Ui(x) 



i e ik fyx ( 2 fc 5 a; 5 / 3 + Qik A x A ' 3 - 6k 3 x -9k^x~- 9i) 9q 2 e~ iJc (k^x - i) 



(26) 



8fc 3 8k 3 
One may easily check that the first order solutions u ~ u + eu\ and v ~ v + EV\ satisfy the Wronskian condition 
up to the order of the calculation. Note that U\(x) vanishes at x = 0, which happens because we dropped all x c 
dependent terms. 

Now, let us check how much the accuracy is improved if we use the first order solutions, compared to the zeroth 
order one. The dominant term in u\ in the high k limit is 0(k 2 ). The first order correction term is 0(ek 2 ) and the 
zeroth order term is Oik). Therefore, the accuracy of the solution is increased by ek, which should be a small number 
for the approximation to be valid. To determine the time when the WKB and the asymptotic approximations will be 
matched, we need to know the size of the error of the zeroth order solution for a given x. The relative ratio of the 
correction term to the zeroth order solution gives the error, 

fkx 4 / 3 . 



E, 



'm{x) 



For high frequency modes, there is a long period where both approximations are applicable: 

k~ 3 < x < e _1 . 



(27) 



(28) 



To minimize the error of the approximation, we choose the time x* and use the WKB wavefunction for x > x* and 
the asymptotic solution for < x < x* . The approximation will be best if we choose the intermediate time x* so that 
the accuracies of the two solution match: 



3 



WKB 



which gives x* = (ek ) . For simplicity, we choose to set x* = 1. Therefore, the small expansion parameter becomes 



e = k 



3/2 



(29) 



Hi/k, which ensures that the present approximation works 



Now, the size of error at x = x* becomes -Easym = k~ 
well for high momentum modes. 

In this subsection, we illustrate the calculation up to O(Hifk) 3 . To do this, we need to get the perturbative solution 
up to the third adiabatic order in e. The second and third order corrections are given by the integrals, 



u 2 (x) 
u 3 (x) 



G(x,y)[u (y)V2(y) + u 1 (y)V 1 (y)]dy, 
G{x 1 y)[u {y)V 3 {y)+u 1 {y)V 2 {y)+u 2 {y)V 1 {y)\dy. 



(30) 



Therefore, up to 0(k ), we get the second order corrections 



u 2 (x) 



32 



-iq 



51q 



7k ^ 



9. + / 9 V 



3q 2 - 1 



16i 



32 



7k 2 x 2 / 3 5k 3 x 



The third order corrections up to 0(k 3 ) are given by the integral, 



u 3 (x) 



Jkx^ 3 13/3 

384 



i (291<? 4 + 96q 2 - 32) / 1 

Q H = — t= + O =tt 

q 21fc^i U 2 



In summary, the solution to order 0(k~ 6 ) becomes 

u(x) ~ Uq{x) + £U\[x) + £ 2 U 2 (x) + £ 3 U 3 (x). 

We have obtained the perturbative solution up to 6 th order in e in Appendix [B] where the nontrivial direction 
dependent correction to the power spectrum appears for the first time. 



(31) 
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C. Matching 



The WKB solution is now matched to the asymptotic solution in the asymptotic region, at x = 1. We illustrate 
the calculation to the order of (Hi/k) 3 . Calculation to higher order to (Hi/k) 6 is given in Appendix [C] Remember 
that the accuracies of the two solutions become equal at x — 1. Therefore, on the whole, the present approximation 
becomes the order by order expansion with respect to ek = k~ 2 = Hi/k. The matching condition is given by 



A+u(l) + A_v(l) = (^^(1) 



A+u'(l) + A.v'(l) = 0Wkb(1) = M 1 + 



3# vW 



2m 



3Hi 



where the phase factor is given by 



The frequency at x* is given by f2» = -| 
respect to f2* and its first derivative are 



i / tidx + iip 



= exp 



(32) 



(33) 



and the ratio of the WKB frequency with 



1 1 3r 2 - 1 1 1 . 

c[x * ] = 1 ~P + F + T^~F + 0( F ) 

rt fx rt f\ f\i 



k V 2k 3 4fc 5 V ; 



(34) 



From the Wronskian relations (fTT|) and (|B8|) , it is straightforward to confirm that the coefficients A + and A_ satisfy 
the normalization condition, to the present order, 



. I 2 I . i2 



4* - (9 X 0)0* 



2iJ 4 fc 3 

As the reference, in the standard isotropic initial universe, the normalized mode function is given by 



iH, 



^isotropic 



/2F 



~ikx 1/3 )e 



l/3\ ifa 1 ' 3 



(35) 



(36) 



D. Power spectrum 



In the limit of the later times x — > 0, using Ui(Q) = with i = 1, 2, 3 • • • , the field behaves as 



<A + + A_) 



3e $» 



2H t k 3 V2 ] V 




WKB 



(!)(«'(!) - «'(!)) -0Wkb(1)K1)-<1))] 



fc u'(l)-u'(l) 



A- 2 



i 1 



in'(x*)\ i /n u(i)-v(i) 

2ft 2 



The explicit value to 0(fc 6 ) becomes 



ffin = 



2H,k 3 



<J>*e 



— /A: 



1 



fc 2fc 2 



|- + i \+i i(l 4 + 12i 1 2 ~ 12 ) ~207g 4 - (96 + 756i)q 2 + 788 

"I - ^ 17^ ^ t . "I - 



2/fc 3 



32k 4 



2016/s 5 



21g 6 + 1206ig 4 - (3780 - 384z)g 2 - 2648z 
8064A: 6 



0I F 



(37) 



Note that the final form of the mode solution gains direction dependent corrections order by order in k. However, 
the power spectrum (fT2|) obtained from the field (f37"l) fails to gain corrections from the standard de Sitter spectrum 
to the present order: 

Hi 2 ^//H,\i 
P =4^ + ((t 
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As shown in the calculation in Appendix [Cj the power spectrum acquires corrections only when we calculate to the 
adiabatic order 0(k~ 12 ). The power spectrum including the corrections becomes 

H 2 j 9(11 - 90r 2 +99r 4 ) 
4^ j 32 

For mode with k ~ lOifi, the relative size of the correction term is of O(10~ 6 ). The isotropy violation at initial 
stage of the universe is not small but its effects on the spectrum for the non-planar modes are suppressed by the long 
duration of inflation and high momentum effect. The correction 0((Hi/k) 6 ) is highly dependent on k to suppress the 
anisotropy effect. 

Before closing this section, it is interesting to compare the power spectrum obtained above with the parametrization 
done by Ackerman, Carroll and Wise Q ■ They expected that the power spectrum in a generic model of inflation with 
broken rotational invariancc will be given by 

P{k) =P (k) {1 +g(k) cos 2 e), 

where Po(k) is the usual scale-invariant contribution, r = cos9 = k^/k, and 9 denotes the angle between the preferred 
direction and the wavenumber vector of fluctuations. Our form of the power spectrum is a generalized form of it. 
There are corrections to the direction independent term and they include higher order terms in cos 9, i.e., cos 9 term, 
not only the quadratic component cos 2 9. Thus higher order terms in cos 9 are not suppressed. 




IV. PLANAR MODES 
A. Power spectrum 

For the planar modes r 2 ~ 0, there appears a region where the WKB approximation may not be valid during a 
period in a > 1. We divide the time into three separate regions divided by the times x\ and cc*. In the region 
x\ > x > x* the WKB approximation is valid. For other two regions, we may find approximate solutions. In the case 
of r = exactly, the adiabaticity parameter diverges in the limit of x — > oo and there is no anisotropic vacuum state. 
The mode r = would behave classically, not quantum mechanically, and will be out of scope in this paper. 

The WKB approximation is valid up to the order of the correction term if 



e(x) :- 



dVt 2 {x) 
dx 



29? {x) 



1 

k \ 



-lex 



1/3 



-lex 



1 + 



tanh ex 



-lex 



r 2 t anh ex 



3/2 



< 1. 



The characteristic behavior of e(x) is shown in Fig. [TJ 




Sitter 



Xtop 



FIG. 1: The behavior of the adiabatic parameter in the case of a planar mode. The left hand side (x oo) corresponds to 
the past. The WKB approximation is temporary violated around x ~ xt op . 
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We are considering the planar mode with r<l. For a small r<e ex , the adiabatic parameter becomes of 

Hi 3 + e' 2ex 



e(x) 



ke~ ex 2 2 / 3 (l — e -2exW3 

We notice, for a large x, that the value of e(x) should satisfy 



173+°^) 



e(x) 



3H, 



< 1 



H 



(39) 



Therefore, e £Xl must be large but should be smaller than k/Hi if the WKB approximation holds. In addition, we 
may also expand e(x) for a large x: 



e(x > 1) 



3£ .l/3 e -2 £ x 



1 



27 - 30r 2 - 5r 4 e 
18 



For a small e _ea: , according to this expansion, the adiabatic parameter increases as e~ 2sx until e~ ex <r<l. It has 
a peak at x = xtop with e~ 6Xtop ~ r when when the 0(e~ 4ex ) term becomes equally important with the dominant 
term. Later it starts to decrease. Therefore, we additionally need 



1 



(40) 



since the point X\ is located after the WKB invalid region. Thus X\ should be large enough so that the large x 
expansion hold, but it should not be too large to satisfy the above two relations. The frequency squared for £i > 1 
becomes 



n 2 = e 2 



-^)+Vl(x), for ear »1. 



where 



2 2/3 fcr 2 2/3 k 4£ 2 q 2 

r = — , q = „ TT , V\ 



3H 



3 e 



4c x 



In this section, we consider the solution only to zeroth order. The corresponding zeroth order solution for Eq. 
is given by the Bessel function: 



a II ■ x,( -\ J -if(<l e for £X > L 

oHi smn(7rr) 

Here we choose the solution so that it becomes an incoming wave form: 

1 



\/2 ■ 2 2 /3fc 3 



:p |zrex + 



for 



X — > 00, 



(41) 



(42) 



1/2 



|) lT . The accuracy of this solution is given by the ratio of the 



where the initial phase becomes e 1 ^ = 

1 1 (1— ir) 

first order solution to the zeroth order one. 

Now we evaluate the accuracy at x = x\ when the WKB approximation becomes valid again. Green's function is 
given by 



Gini(x,y) 



3H 



is 



(43) 



Noting that qe £Xl 3> 1 because of Eq. (|3"9")l . the error can be estimated to be 







00 / 


4>o 







G- m i(x,y)V 1 (y)^o(y)dy 



Sexi 



H, 



Let us determine the cutoff time x\ so that both of the solution (|4T|) and the WKB solution arc equally valid. The 
WKB solution is valid if e < 1. For ei> 1, the accuracy of the WKB approximation is 



E 



WKB 



( 3 H 



V2 2 / 3 k (r 2 + e - 2ex fl^ 



J2ex 
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where we use Eq. (j40|) . The WKB approximation is valid during 

1 < e sx < min(-^-,r -1 ). 

Both approximations are valid for this region. Setting the two errors be of the same order, we choose x = x\ to be 



=Cl\r=r 



Hi 



3/5 



where c\ is an O(l) constant. It shows that the time x\ is dependent on momentum but independent of directions. 
Note that this value e £X1 should be smaller than r _1 because of (|40|) . Therefore, we have constraint for r: 

r « (ft/fc) 3/5 - 

The accuracy at x = x\ is given by 0(ft/fc) 4 / 5 . 

During x* < x < X\, we use the WKB approximation where is given in Eq. (|29[) . The WKB solution is 



ft 



exp • 



|«y dx'Q(x')} + -^L=cxp{ -ij dx'n(x')} 



At x = xi, to the zeroth order in e, the matching condition becomes 



where 



fl 1 



f> 2 ~e 2 (f 2 



(f>o(ti) - i 



v2sx 



■)=e' 



2 4 / 3 fc 2 
9Hf 



^o(ii) + i 



mi) 



(44) 



(45) 



1 /ft\ 6/5 \ 2 4 /3 /HA n/5 



The argument of the Bessel function satisfies ge ea:i 
asymptotic form for the Bessel function: 



til ( Jl 

3 I Hi 



k 

2/5 



9 C 2 



^> 1. Therefore, at x = x±, we may employ the 



lim Ji a (a;) 

a;— »oo 



cos a; 



7T ICLTT 



and find 



The coefficients B+ and _B_ become 



2 . 7T ir7r, 

cost oc " — — -I- ) 

6ft sinh(Trf) V irqe- eXl v 4 2 ; ' 



(46) 



ft 



e e 



:(f-fe-"i) 



3ft Vl - e" 2 ^ 



£ e e H -i 



3ft Vl — e 



(47) 



Note that (ft/fc) 2 < f < (fc/ft) 2 / 5 and ft satisfy |ft| 2 - |ft | 2 = e/(3ft). 

For x -C a;*, the WKB approximation is not valid. However, for x < x* the asymptotic solution (|37[) works. For 
< x < x*, the isotropic solution is valid and 

4>o = A + u (x) + A_v (x) . (48) 

Matching the asymptotic solution with the WKB solution similarly to the previous section, we obtain 

ft?\* 



a; -)• 



o>~(^ + ^-) = (§r)V* 



— 7rr i^IM 

e e 



where ^(/c) = fcx* — J** Q(x)dx + qe 5X1 — j. This leads to the power spectrum of the form: 



'-(§)'( 



coth 7rr 



cos (2*) 
sinh nf 



(49) 



(50) 
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The explicit value of the correction term becomes order of 10 -3 for f ~ 3. 

At a glance, the phase factor VP looks sensitive to the choice of the matching points x* and X\ explicitly. However, 
it can be confirmed that it is independent of them. Now, we compare two sets of the matching points (x±,x*) and 
(xi,x*), and take their difference to obtain 

* - * = fc(a-y 3 - ic|) + q(e' exi - e~ eSl ) - f ' Vl{x)dx + f * Q(x)dx 

~ k(xl /3 ~xl) + q(e' EXl - e~ ££l ) - f ' -^dx + f ' qee^dx = 0, 

where we have used the values of the frequency squared around x — x* and x = x\. 

After averaging the power spectrum (|50[) over the plane, i.e., integrating over the momenta k\ and &2, the second 
term of the spectrum vanishes and the main contribution comes from the first term, which is approximated as 

(f) 2 coth M = (|) 2 -th(^cos.), (51) 

where we have used the definition of f and r = cos 9. Thus in the planar limit, the deviation of the power spectrum 
from the ansatz in Ref. [jfl is quite clear. It makes it impossible to compare observations through the analysis based 
on the ansatz in 0, 120U22] . 

B. Physical bounds for r 

Our analysis in this paper indicates that the power of the scalar field fluctuation has a negligible difference from 
the standard inflation in the non-planar directions, but it has a sharp peak around the plane 9 = w/2. In other words, 
the power is strongly localized around the plane 9 = tt/2. While such a feature is much different from the often used 
phenomenological ansatz for the anisotropic power spectrum Refs. 0, H3|, it also implies a band- like structure in the 
CMB temperature map, around which the averaged temperature is much different from that in the other direction. 
Recently, it was reported that giant ring structures which have an anomalous mean temperature profile were observed 
in the CMB sky (23l . |2~H . Thus it would be very interesting to compare the predictions in our model with such a new 
observational feature. 

In the rest of this section, we discuss the possible bound on the parameter. In particular, while the power spectrum 
in the planar mode Eq. (fBT)]) diverges in the r = limit, we also find the purely planar mode r = cannot be 
quantized. It implies that we need to introduce a cut-off on the r to be consistent with theoretical and/or observational 
requirements. From the above analysis, we have found that for a large f > 2, the power spectrum (|50|) is not much 
different from that of the de Sitter (the difference is too small; the difference is 0(1O~ 5 ) for r = 2). Thus we focus on 
the case f < 2. 

A bound is obtained from the validity of our method. The present approximation holds for (k/Hi) 2 / 5 ^> f ^ Hf /k 2 . 
The amplification will be maximized for a small f. By defintion of f, for a given fc, it is equivalent to put the bound 
on r as (k/Hi)~ 3,>5 ^> r Hf/k 3 . This bound is from the validity of our approximation. But it can only give a 
/c-dependent bound on r. However, since fc and r are independent variable in the Fourier space, a fc-independent 
bound may be obtained from the consistency with requirements of theory and observation. Now, we would like to 
check whether our scalar field solution naturally satisfies the required bound or not. 

From the theory side, we may check the backreaction of the scalar field may modify the background evolution in 
the infinite past where the scalar field is quantized in the anisotropic vacuum. Since the UV divergence may be cured 
by the standard renormalization, here we focus on the IR behavior (For UV, here we simply put the UV cut-off k uv , 
which may be chosen close to the Planck scale m p ). Now we check the sensitivity of the scalar field backreaction to 
the IR cut-off ki r for /c-integration and the fiducial cut-off < r$ -c 1 in the r-integration. If the energy density is 
very sensitive, it would lead a bound on ki r and tq. In the Appendix [D] we explicitly compute the energy density of 
the scalar field p after introducing the ultra-violet and infra-red cutoff and the leading order result in the infinite past 
is given by 

where X = e~~ 2x . As seen here, there is no infra-red divergence, leaving no bound on the IR part and one may 
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naturally take tq — } and ki r — > limit. 1 

Finally, we turn to the observational side. One of the most important requirements from observations is that the 
anisotropic contribution on the power spectrum should not exceed the isotropic one. In the Appendix [Ej we have 
computed the corrections in the CMB temperature anisotropy. Note that although in this paper we have not discussed 
the metric perturbations, the previous works show that our scalar field perturbations share the same features with 
the metric perturbations (l2l . [l5[ and here we proceed with the power spectrum obtained for the scalar field. After 
some computations, for a given £, the ratio of the corrections from the anisotropy AC« m to the standard isotropic 
angular multipolc moment is given by 



AC, 



ttm 



Cr 



C(£,m) [(kirf]*) lnr 



(52) 



where ?y* is the comoving distance from the last scattering surface to the present time and C (£, m) is given by 
Eq. (|E29|) . Note that the above result is valid when £ + m is even. Avoiding the large anisotropic contribution, 
ACu m /Ci < 1, gives a lower bound on ro for each £ and m: 



r > exp 



1 



C(l,m){k ir r*) 



(53) 



Among the bounds obtained for various £, we should choose the strongest bound. For m = 0, for example, C(2, 0) ~ 
0.267, (7(10,0) ~ 0.0576, (7(20,0) ~ 0.0284, and hence for a larger even £, C(£,0) becomes smaller. Thus, we should 
impose a bound on ro from one of the lowest £: focusing on £ = 2, 



rp > exp 



3.70 



(54) 



If we put fci r ?7„ = 1, the right-hand side of Eq. (|54[) becomes 0.0246, which is reasonably small. As shown in the 
Appendix. 02 Hit)* > 2. Choosing the IR cut-off ki r < Hi, the above lower bound can be naturally obtained. In 
other words, the lower bound can be made less serious because we may choose the parameter ki r freely. If the giant 
ring or any other specific effect would be observed in the future, it would finally determine the value of Tq. 



V. CONCLUSION 



In this article, we have reinvestigated the quantization of a masslcss and minimally coupled scalar field as a way 
to probe the signature of pre-inflationary background anisotropy in the spectrum of cosmological perturbations. We 
have considered the gravitational theory composed of the Einstcin-Hilbert term and a positive cosmological constant. 
For simplicity, we have assumed that the late-time inflationary stage is exactly described by the de Sitter solution 
and ignored the dynamics of the inflaton field, which can be justified by the slow-roll approximations as shown in 
Ref. [l2j. The initial geometry we have considered is the axially symmetric Bianchi I metric which evolves to the 
de Sitter spacetimc due to the positive cosmological constant. The motivation to consider the initially anisotropic 
Universe is two-fold: The first is that even if the current Universe is almost isotropic, it does not mean that the 
Universe is isotropic from the beginning. In fact, Wald's no-hair conjecture ensures that in the presence of a positive 
cosmological constant an initially anisotropic Universe exponentially approaches a de Sitter spacetime at the later 
time under the strong or dominant energy condition. Therefore, it would be more generic that the initial Universe 
is anisotropic and one expects that remaining initial anisotropy is imprinted in cosmic observables. The second 
motivation comes from the recent observations by WMAP satellite which almost confirmed the predictions from 
inflation, during which the Gaussian and statistically isotropic primordial fluctuations arc produced. But, several 
groups have reported the possibility to detect the remnants of the preinflationary anisotropy through the fluctuation 
spectrum (See [25[ and references therein). 

We first dealt the non-planar modes. We have divided the time into two periods, before and after the mode- 
dependent time r» ~ jj- log We used the WKB approximation for the first part (t < t„) and used the asymptotic 
approximation based on the de Sitter solution for the next part (r > r*). At t = t», the two approximations have 
the same accuracy with an error of order Hi/k. We illustrated the calculation up to the order of (Hi/k) 3 in Sec. IIII CI 



1 For a different point of view, see [l7l | ■ 
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and up to the order of (Hi/k) 6 in Appendices A to C. We have shown that, at each order, the asymptotic value of 
the scalar field acquires corrections from that of the de Sitter form and has the form, 

<t>fin* = (A* + B k i)e- i VV^, 

where and B k are real constants which are dependent on the momentum. Since the phase part is written as a 
multiplicative factor, |</>fi n | 2 = A^ + B^ gains no momentum dependent oscillatory behavior. It was also shown that 
the power spectrum of the scalar field acquires non-vanishing corrections only when we execute the approximation 
up to 6 th order. Hence, the direction dependence appears only at the order 0((Hi/k) 6 ). If we want to observe the 
direction dependence from the power spectrum (other than the planar mode) , the wavelength of the fluctuation due 
to the modes with Hi < k CT it < k < lOHi should be located in the observable range of the universe (1 — 10 4 )Mpc. For 
modes with k > lOHi, we can not detect the effect with the present equipment because it is too small (the fluctuation 
will be smaller than < 10 -6 ). For modes with k < Hi, the approximation we adopted in this paper is unapplicablc 
and we cannot apply the result In addition, it also restrict k CIlt < 10i2j, which gives TV < 64 + log e (10) ~ 66 in 

the case of GUT scale inflation. In summary, if the power spectrum is observed in the previous form, it restricts the 
number of e-folds to the narrow range 64 < N < 66. If the number of e-folds is larger than 66, we may claim that 
the effect of anisotropy will be unobservable within the present accuracies of observations. If the number of e-folds 
is smaller than 64, some of long wavelength fluctuations, which are of the order of the scale of our universe, cannot 
be dealt with the present adiabatic approximation. The two-point function in this regime was numerically analyzed 

For the planar mode, we have obtained essentially the same result as that in our previous analysis |12Lll5|. For such 
a mode, the temporal breaking of the WKB approximation relatively enhances the effects of the primordial anisotropy 
in the power spectrum. For j^- ~ 3, the deviation from the ordinary case becomes of order 10 -3 . We also discussed 
the observability of the planar mode contributions. It was shown that there exists a sharp peak around the symmetric 
plane. The giant ring structure in the CMB sky map can be related with this result. Note that our result is applicable 
even for a small &3 ~ Hi/10. In this case, the power spectrum increases by a factor 6 to the standard de Sitter case, 
which was not observed. There are three possibilities to cure this defect. The first is that there may exist a lower 
bound on the momentum orthogonal to the plane, which was discussed in Sec IV B. The second possibility is that 
the large energy density of the scalar field reacts to change the background geometry. The third possibility is that 
the gravity theory can be modified at the ultra high energy scale so that the anisotropy may not grow indefinitely at 
initial times as in the case of the Eddington-Inspired gravity (26j . 

The properties of the non-planar mode function also suggest that it would be very difficult to obtain distinguishable 
amounts of the primordial non-Gaussianities in the present model. The bispectrum is very sensitive to physical process 
before, during and after inflation. Assuming three comoving momenta (i = 1,2,3) which forms a triangle in the 
Fourier space, • kj = 0, interactions just before and after the horizon crossing (i.e., classical interactions) give rise 
to the equilateral-type (ki ~ k% — &3 where ki :— |k»|) and local-type (ki <C min(/c2, ^3) and permutations) bispectra, 
respectively. On the other hand, the existence of an excited state in the past contributes to the folded-type bispectrum 
(ki + &2 ~ £>3, and permutations) Since after the isotropization the evolution is the same as the conventional 
single-field slow-roll inflation, only the possible origin of the distinguishable non-Gaussianities may be the folded-type 
one due to the modification to the anisotropic vacuum. The absence of the negative frequency mode, however, implies 
that a negligible difference from the standard case for the non-planar high-momentum modes. 

For the planar mode, since the contribution of the negative frequency mode is less suppressed than for the non-planar 
mode, the amplitude of the folded-type bispectrum may be nontrivially different from that of the standard case and 
have features of modulations. In the present model, however, the overall amount of the primordial non-Gaussianity 
is already suppressed by the slow-roll parameters, and hence it would be still very difficult to detect their difference 
with the precision of the near future observations such as the Planck satellite. Non-Gaussianities may be enhanced 
in the inflationary models with noncanonical kinetic terms for the inflaton fields. It would be interesting to observe 
the process of the isotropization in such a generalized inflation model. 

Acknowledgement 

HCK was supported in part by the Korea Science and Engineering Foundation (KOSEF) grant funded by the Korea 
government (MEST) (No.2010-0011308) and the APCTP Topical Research Program (2012-T-01). MM is supported 
by Yukawa fellowship and by Grant-in- Aid for Young Scientists (B) of JSPS Research, under Contract No. 24740162. 
MM is also grateful to the School of Liberal Arts and Sciences, Korea National University of Transportation, for the 
hospitality. 



14 



Appendix A: Higher order WKB approximation at early times 

We need to solve Eq. (fT6|) iteratively. To zeroth order, the frequency squared is given by 

r> 2 = n 2 . 

Including the first order correction, it is given by 



1 - 









= n 2 


\nn J 





1 - 



2ft 3 4 V Q 2 



which is correct to the order of e n+1 . Note that the first order correction terms are of order 0(e). In this way, the 
higher order corrected frequency squared is given by the iteration, 



6 2 — O 2 



1 - 



3 / a 



2n 2 n n 4 I na n 



To have accuracy up to 0(e 5 ), we need to get ft 2 iteratively. For ex <C 1, through the successive iterations, we 
obtain 



1 1 3r 2 - 1 1 \-r 2 549r 4 - 498r 2 + 133 3r 2 - 1 -45r 4 + 41r 2 
!__ + _.' 



x = l 



k 2 fc 4 6fc 5 k 6 k 7 72k 8 
-14175r 6 + 19629r 4 - 7191r 2 + 485 -549r 4 + 486r 2 + 35 



2/s 9 



3/fc 10 



324/c 1 



96k 1 



(Al) 



In addition, we obtain 
Q'(x 



bq 2 9q 2 6<? 4 + ^--f 3 (45g 4 + 12g 2 - 4) 11 (405g 6 + 162g 4 - 54g 2 - 
1+ 2F + 4F + p + 324fc 9 



5 (2079g 6 + 756<? 4 - 252g 2 - 32) / 1 
48fc n + °{W 2 



(A2) 



Appendix B: Perturbative calculations in the later times 

We have used the iterative perturbation method. To solve Eq. (fT8|) . we first set the solution in series form as 

oo oo 

u = ^2e n u n (x), v = e n v n (x), 



n=0 



n=0 



where u and v are two independent solutions of the equation of motion (|18|) and are complex conjugates to each other. 

The zeroth order solution is given in Eq. (f2"2"|) . Once we obtain the solution up to n th order, (n + l) th order correction 
can be found, by solving the following differential equation, 



dx 2 9a; 4 / 3 



i = - u m(y)v n - m (y), 



(Bl) 



with the source term becomes 

1 Q-r.l/3' 



Vo = 

2fc 2 (r 2 - 7) x 11 / 3 
10935 ; 



2fc 2 r 2 x 2 / 3 



Vo 



8k 2 x 5 / 3 



27 729 

2fc 2 (I89r 2 + 170) x 14 / 3 
6 = 2066715 



Vi = - 



k 2 (ll + 30r 2 )a; 8 / 3 



10935 



(B2) 
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Therefore, the (n + l) th order solution can be obtained, by using the Green's function (|25j) . as 



u n +i{x) = dyG(x,y) 



^2 u m(y)Vn-m(y) 



By using this method, we get 



it 2 (x) 



+ e 



-ife ^/x 



9iq 4 a; 5 / 3 27 (g 4 x 4 / 3 ) 27ig 4 x 9z (99g 4 + 24g 2 - 8) 9 (99g 4 + 24q 2 - 8) 



32k 



32k 2 



32k 3 



64k 5 



64k 6 



1 , (459<7 4 + 96q 2 - 32) x s / 3 k 2 1 , . „ . 7 ,,- 

i<7 4 x 3 fc 3 + ^ ^ ^ + — « (99g 4 + 24q 2 - 8) x 7/3 fc 



32 



2016 



+ 



— (99<? 4 + 24 g 2 -8) x 2 
40 v 1 H ' 80k 

9i (99q 4 + 24q 2 - 8) $x 9 {99q 4 + 24q 2 - 



9i (47g 4 + 12g 2 - 4) x 5 / 3 3 (45<? 4 + 12g 2 - 4) x 4 / 3 3i (45q 4 + 12q 2 - 4) 



16k 2 



16k 3 



64fc 5 



64fc 6 



(B3) 



w 3 (x) 



^— ik ^fx 



9<? 6 x 3 i<z 2 (459q 4 + 96g 2 - 32) x 8 / 3 9 (q 2 (99q 4 + 24<j 2 



.7/31 



9iq 2 (99g 4 + 24q 2 



256 



+ e 



— (Vx 13 / 3 ) fc 4 
384 V J 



1792fc 896fc 2 320fc 3 

i<7 2 (291q 4 + 96<? 2 - 32) x 4 fc 3 (102789<? 6 + 40176g 4 - 13392? 2 - 1792) x n / 3 fc 2 



8064 



362880 



i (50193g 6 + 19224g 4 - 6408g 2 - 896) x w ' 3 k (-273807*? 6 - 101952q 4 + 33984g 2 



30240 



34560 



i (284607g 6 + 104352 g 4 - 34784 g 2 - 4480) 
8960fc 

i (288387(7 6 + 105192? 4 - 35064 9 2 - 4480) x 2 



n 8/3 



(288387g 6 + 105192^ - 35064g 2 - 4480) 
2688fc 2 



,7/3 



960fc 3 
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(B4) 



-ik %fx 



3^ g 8 x 13 / 3 fc q 4 (291g 4 + 96g 2 - 32) . 



ik <yx 



l x 17 / 3 k 5 (q 4 (477q 4 + 192g 2 - 64) x 16 / 3 ) k 4 



6144 



1024 7168 

i (484623(? 8 + 247752 9 6 - 71064g 4 - 20224g 2 + 1280) x 5 fc 3 
10160640 

(30240459g 8 + 15875352g 6 - 4410504 9 4 - 1553408g 2 + 141824) x 14 / 3 P 

66044160 

i (5658255g 8 + 2920077? 6 - 812511q 4 - 278256g 2 + 25712) x 13 ' 3 k 

1572480 

(64266939g 8 + 32444658g 6 - 9097542 g 4 - 2976096g 2 + 269216) . 

2661120 



129024 



(B5) 



u 5 (x) 



q 10 x 7 k 6 + iq 6 (2133 9 4 + 960q 2 - 320) x 20 / 3 P 



ik \/x 

' 122880 ' 7741440 
(q 2 (853821q 8 + 519264g 6 - 127008g 4 - 55808g 2 + 5120) x w ' 3 ) k 4 

162570240 

s (175688271 9 10 + 114725160? 8 - 25188408g 6 - 16032384g 4 + 1286400 9 2 + 186368) x 6 P 



u 6 (x) 



ik tyx 



2377589760 

^12^,25/3^,7 q 8 ( 993(7 4 + 48Qg 2 _ 160 ) ^8^6 
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(B6) 
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The independent solution v(x) is nothing but the complex conjugate of u(x). The Wronskian can also be checked to 
give 

2ik 3 



+ o(k- w ). 



(B8) 



Appendix C: The asymptotic form of the scalar field 

In the limit of the later times x — > 0, the wavefunction behaves as 



"liu 



^ Q = -(A + +A_) = |^wkb(1)K(1) -«'(!)) -^Wkb(1)(«(1)-«(1))] 



3^jM / k u'jl) -v'{l) I ijVfa) \ Kl u(l)-v(l) 
2H t k 3 ^2 I V ^* fc 2 V + 2^2 ) Y k k 



By using the perturbative results in the previous Appendices [A] and [B] we get 

% (q i + \2iq 2 - 12) 



1 _i_ l£l 

2 4 



2# ! /fc 3 e ifc _ . 1 
^— ¥ -</>fi„-*+ I + -^ 2 p 

-207<7 4 - (96 + 756t)g 2 + 788 21q 6 + 1206ig 4 - (3780 - 384i)< ? 2 - 2648i 



2016/c 5 8064fc 6 
i (825g 6 + (480 + 74142i)g 4 - (19060 - 18432i)g 2 - 18744i) 

40320F 

i (945g 8 + 446904^ 6 - (11016648 - 262656i)g 4 - (2654208 + 3262752i)q 2 + (2472336 - 28672i)) 

5806080P 

13905g 8 + (8640 + 3002616i)<7 6 - (37893384 - 926208i)<? 4 - (9103104 + 3483936i) 9 2 + (4621968 - 28672i) 



1 



5806080fc 9 

1323 9 10 - 2928798ig 8 + (831208392 - 2149632i)g 6 + (338204160 + 1072579536i)g 4 



162570240fc 10 

(212753520 - 255210496i)g 2 - (17461248 + 124990304i) 

i 



1 410319(7 10 + (262080 + 260683110i)<7 8 - (14181596520 - 142566912i)g 6 

21134X3X20& 

-(5225036544- 163571184i)q 4 + (3041923248 + 133211136i)g 2 + (226996224 + 471614432i; 



5733g 12 + 41102100ig lu - (23757489540 - 34085376i)<f - (11283840000 + 1401723182496i)g b 



16907304960fc 12 

- (232347906000 + 516373012480i)<? 4 - (57447653376 - 183565086016i)<? 2 + (15556834880 + 22417833984i) 
+0(k- 13 ). 

Notice that this is of the form </>fi n k = (Ay L +B] i i)e~' lk where and B^ are real numbers. Therefore, |0fi n k| 2 = A^+B^ 
has no oscillatory behavior with k. 



Appendix D: Quantum backreaction in the infinite past 



In this Appendix, we discuss the quantum backreaction of the scalar field in the infinite past. We rewrite the metric 
Eq. ([5]) into that written in the new coordinate x, defined by sinh(3iJir) = — ij- - . In the x- frame, the metric can be 
rewritten as 



ds 2 



1 \ 2 



where 



o**(J_ 



o a ( x ) 



dx 2 



e~ 2p {dx{ + dxi) + e^dxl 



1 



sinh 1/3 ; 



e /3(*) =e -f*. 



(Dl) 



(D2) 
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We are interested in the early time behavior of the backreaction and adopting the leading order WKB approximation, 

ji / cte'fij, 4>k,x — ifl<t>k- (D3) 



1 1 



Decomposing the scalar field into the Fourier modes as Eq. © , the energy density of the scalar field is given by 

1 



P 



2(2tt) 3 



d A k 



— 6a / 



-(3H i ) 2 |^r + e- 2 «(e 2 ^i + e- 4 ^l)|^| 2 ' 



(D4) 



Substituting a(x), /3(x) and mode functions, the energy density is rewritten in terms of the integration of the Fourier 
mode 



p{x) 



1 



d*k 



2(2tt) 3 

Substituting Eq. ((D3)) into Eq. (|D5|) we get 



(3H t ) 2 sinh 2 x\<p Kx \ 2 + smh 2 l z x[e- 2x ' z kl + e 4 */ 3 fc 2 ) \cf> k \ 



(D5) 



,,2.r 



4(2tt) 3 \3H,, 



d A k- 



(3^) 2 (|) (1 - e- 2x ) 2 n 2 + k 2 (l) 1/3 {l - e- 2 *) m (e- 2 *{l - r 2 ) + r 2 )] . (D6) 



For the further convenience, introducing a new coordinate X := e 2x , and using d 3 k = 4irk 2 dkdr (here, we have taken 
the symmetry between the north and south hemispheres into account) we write the energy density as 



P{X) 



TlX- 1 / 1 



dkk z I drU k (X,r), 



where we explicitly introduced the ultra-violet and infra-red cutoffs explicitly. Here, 



U k (X,r) := I 



(1 - X) 2 n 2 + k 2 (\) 1/3 {l - X) 2/3 {(1 - r 2 )X + r 2 ) 



(D7) 



(D8) 



Now the UV divergence is cured by the UV cut-off k uv which may be chosen close to the Planck scale m p . In the rest 
of this Appendix, we focus on the IR divergence. 
Following the definition 



n 2 (x) 



k \ 2 2t [(l-r 2 )X + r 2 ] 

m) (i-x)4/3 



in the infinite past, the energy density is given by 

p(X -> 0) 



3 



-tttt- , Q / dkk^ I drr. 
2V3(27r)3 J kir J ro 



(D9) 



(D10) 



Hence there is no IR divergence both in the k- and r- integrations. 



Appendix E: CMB anisotropy 



In this Appendix, we evaluate the CMB temperature anisotropy 4jr in our model. The power spectrum of the 
primordial fluctuations is in general given by 

<(J(k)<J*(k')> = (27r) 3 P(k)(5 (3) (k-k'). (El) 

Here, we assume that the power spectrum of the primordial perturbations including the effects of the primordial 
anisotropy is given by 

P(k)=P o (A:)(l+0(fc ) O fc )), (E2) 
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where fl k denotes the solid angle in the Fourier space and Po(k) is the power spectrum in the isotropic case. Note 
that we have not worked on the metric perturbations, but following the results obtained in (l2l . [l5j we can expect 
that the form of Q is not modified much from the case of the scalar field even in cases of the metric perturbations. 
The CMB temperature anisotropy can be expanded in terms of the angular multipole moments given by 

dnY em {Q)^(n), (E3) 

where Y^ m (0) is the spherical harmonic function satisfying the normalization condition 

dQY em (n)Y e , m ,(Q) = 5 u >5 mm >. (E4) 
The expectation value of the angular multipole moment is given by 

{a, m a% m ,) = j dQdn'^(n) S ^(n'))Y; m (n)Y eim ,(n'). (E5) 

The primordial perturbations are related to the observed CMB anisotropy by 

= / (^E(^)H) L ft(^^( k )^(fc), (E6) 

where 0l(&) is the transfer function and <5(k) is the primordial perturbations in the Fourier space. Here, we ignore 
the Integrated Sachs- Wolfe (ISW) effects and focus on the Sachs- Wolfe (SW) one, the transfer function is given by 

47T 

&t(k) « -yj<(fo?.), (E7) 

where 77* := 770 — T)i s , Vo and 771s are the conformal times of the present universe and of the photon last scattering, 
respectively. 

The angular correlation function of the CMB fluctuations is given by 



(f(fi)fV)) = / (gs E 1} {2L 2 l) (-i) L - L 'P L (k ■ n)P L 4 ■ n>)PQL)e L (k)Q L ,(k). (E8) 

Using the relation of the Legendrc function and the spherical harmonics 

PL{k ' i)= 2LTl E Y£ M (Sl)Y LM (Cl k ), (E9) 



A/=-L 

the angular correlation function is given by 

d 3 k 

L,L',M,M 



(^(n)^(fi')} = / 7^ E {-i) L - L 'Y L M{mLi^ k )Yt IA ^n')Y L , M ,{n k )p{^ 

"* V ' T. T.i M Ml 

Substituting (|ET0|) into (|E5|) . 



dkk 2 

a lm a\, m ,) = {-if-"' I — - I dn k Y; m (n k )Y e , m i(n k )p(k)Q e (k)e e ,(k), (eh) 



(2tt) 3 

where we have used the normalization condition Eq. (|K4|) and d 3 k = k 2 dkdft k . Using the power spectrum (|E2[) , 
(at m a* e , m ,) = S eel 5 mml ^-(4ny 2 J dkk 2 P (k)Q e (k) 2 

+ (-*)<-<' {^Y 2 \ J dkk 2 p (k)e e (k)Q e ,(k)( I dn k Y* n (n k )Y e , m ,(n k )g(k,n k ))- (E12) 

Since the correction is negligible for the non-planar mode, we may focus on the planar mode. It is adequate to 
approximate the whole power spectrum with that for the planar mode. Since cos(2 , I') term in the power spectrum 
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Eq. (|50p is oscillating, we may ignore the angular power spectrum because of the rapid oscillation. Thus the function 
Q is given 

0(fc,n fc ) = coth(^^|cos0 fc |) -1. (E13) 

Considering the SW approximation, Eq. (|E12|) becomes 

2 



(a lm a\, m ,) « — ^Su'Smm' j dkk 2 P (k)j^(kri„) 2 

+ {-if- 1 ' j dkk 2 P {k)jt{kr]*)je(kri*) f dn k Y; m (n k )Y e - m ,(Cl k )f(k, |/i fe |)} , (E14) 



where fi k = cos6 k and f(k, \fi k \) := coth ( ^gjj^ \^ k \j - 1. 
Noting the definition of the spherical harmonic function 



F,„(S! t ) = (-irW^- 
the integration along the azimuthal (^-direction reads 

dn k Y; m (n k )Y e , m ,(n k )f(k,\n k \) 

Om.', 



(E16) 



We then perform the integration along the /i^-direction. Noting the parity transformation of the Legendre functions 
we can rewrite 

/ dfM k IT(ti k )PF(n k )f(k, = fl + (-l)^'l / drPr(r)P?(r)f(k,r), (E17) 
J-i L J Jo 

where r := %■ = \^ k \. Substituting into (|E14|) 



M, m .) « ^{fa. / dfcfc 2 P (^(^) 2 + [1 + (~1) W '] ^ (2< + l)(2f + 1) | - gj| - g 

x / dkk 2 P (k)j e (kri,)j t (kv*) f drP™(r)Pp(r)f(k,r)}, (E18) 



where the lower bound on the r- integration < ro <C 1 is added as a cut-off to avoid the divergence from f(ro — > 0, fc). 
The nonvanishing averaging is obtained for m = ml . 

From now on, we assume the exactly scale-invariant power spectrum 

k 3 P (k)=A s , (E19) 

where A s is the amplitude of the primordial perturbations. For later convenience, we decompose the averaged multi- 
pole moments into the isotropic part and the anisotropic corrections as 

(ai m a* e , m ,) = 5 mm 'iCe.Sw +AC«/ OT (ro)|, (E20) 

where the isotropic part Ce is given by 

Ce = ^f dkk 2 P (k)j e (k^) 2 , (E21) 

and 



97T y (£ + m)!(r + m)\ 

dkk 2 P (k)j e (k v ,)j e ,(h h ) / drPP(r)PP(r)f(k,r), (E22) 

ir ^ r O 
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respectively. Here ki r is the IR cut-off for the anisotropic mode. Using (|E19I) . the isotropic part reduces to the 
standard result 



1(1+1) 
2tt 



C, 



18tt 2 ' 



(E23) 



which gives the SW plateau seen for £ < 100. 
For a given £, we take the ratio 



dy 



Mw.H i )j i ,(yri.H i ) f drPr(r)PP(r) ( coth ( 



2 2/3 



7r?/r 



(E24) 



where we have changed the integration variable y = -i-. An observation is that there is no contribution from the 
anisotropic part if I + £' is odd. Thus we focus on the case £ + £' is even: 



1) £ = even and 



even 



2) £ = odd and f = odd. 

In the integrand in Eq. (|E24|) . the dominant contribution comes from the region of r ~ 0. In case of 1), if to = even, 
Pp(r) ~ const for a small r, while if m = odd, P™(r) ~ r for a small r. In case of 2), if to = odd, P™(r) ~ const for 
a small r, while if to = even, P™(r) ~ r for a small r. Thus the integrand of r is divergent as r _1 in the cases 1) with 
an even m and 2) with an odd to, while it is proportional to r in the cases 1) with an odd to and 2) with an even to. 
Thus we may focus on the first case, to avoid too large anisotropic contribution to the CMB angular spectrum. 

We now focus on the case of £' = £. From Eq. (|E24j) . 



AC« m (r ) _ v ( 2 £ + l)(£-m)\ 



°° dh r 1 / 

YMfoh) J drPPir) 2 (coth 



2 2/3 



ixyr 



Ci K ' (£ + m)\ 

Thus when £ + m is even, the leading order contribution from the r-integration is given by 



1 



(E25) 



j™ h ^(H^yf jy r Pr(rr (coth (^^) - l) 



3 x 2 2m lnr ( 7 1 



2 2 /3r(I(l - £ - m))T(l + - m)) 



(E26) 



where 



-^ji(kirv*yY 



- 4 



(2* + 3)(2£+l)(2*-l)r(§+*) 

1 . „i r„ 1 „ 3 



{(kr»7.)r(|+i) 5 



" 1 -'(fc 4 , ? 7*) 2£ (3 + 8£ + 4£ 2 ) 2 P 3 [{* - 2' 1 + + 2 ' £ + 2' 2 * + 2 >' -( fc ^*) 2 ] }' ( E27 ) 



Note that for £ > the second term in Eq. (|E27[) is negligible. Hence from Eq. (|E25|) . the leading order ratio 

becomes 



ACt 



C(£,m) [(fc ir 77*)lnr x ] , 



where 



C(£,m) 



£(£ + !)(£ -m)\ 2 1 / 3 tt 



3x2 



2 in 



(£ + to)! (2£ + 3)(2£ - 1) r (i (1 _ e _ TO )) 2 r (l + \{i - m ))' 

To avoid a too large anisotropic contribution, AC« m /CV, gives a lower bound on ro 

1 

?'o > cxp 



C(£, m)(k ir r)* 



(E28) 



(E29) 



(E30) 
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Before closing this Appendix, we discuss the possible value of 77* . For simplicity, we ignore the contribution of the 
dark energy and assume that the universe is matter-dominated after the last scattering. Representing the proper time 
in the matter-dominated universe by r m , the scale factor is given by 



a m (T m ) = a 



T m.O 



(E31) 



where r rn ^ is the present time and ao is the present size of the universe. The conformal time, defined by dr] = 
), is now given by 



V = 



o 2/3 1/3 
6 T m,0 Tm 



The comoving distance to the last scattering surface is given by 



Vo - m s 



3r m 



1 



7"m,ls 
T m ,0 



1/3 



3r TOl o 
a 



(E32) 



(E33) 



where 



k < 1. 



The modes which contribute to the SW effects reenter the horizon after the last scattering r rn j s < r m < r m ^ - 
We assume that the same modes cross the horizon during inflation at the time r, where the scale factor of the 
universe is given by a; n f(r) = e * T , and the horizon exiting time during inflation r and the re-entry time during the 
matter-dominated era r m is related by 



where H„ 



k = a- mi {T)Hi = a m (T m )H m {T m ), 
a' m (T m ) / a m (T m ) = 2/(3r m ). Since a m (T m )H m (r n ^ 



(E34) 



2/3 1/3 ? 



3tJ „t. 



we obtain 



Hi 



O'jriHfj 
flint 



2a 



-Hi 



n 2/3 1/3 
6T m,0 Tm 



(E35) 



which relates r to r m . Assuming the mode which reenters the horizon today exits the horizon at r = tq in the early 
universe, by definition 



Using (|E33[) . we obtain 



H 



2a 



2a 



3r, 



e -ff^ 

m.O 



-HiT 



3r, 



m,0 



3t to ,o 
a 



2e~ H ' T0 . 



(E36) 



(E37) 



Note that the dependence on ao is canceled out. If tq = 0, i?;??* = 2. Assuming the largest mode we observe today 
crosses the horizon just at the onset of inflation, it may be natural to choose To ~ 0, leading to if;?/* ~ 2. Thus the 
value of -ffi?/* should be taken around 2. 
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